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Abstract 

The number of independent sets is equivalent to the partition function of the hard-core lattice 
gas model with nearest-neighbor exclusion and unit activity. We study the number of independent 
sets mdfi{n) on the generalized Sierpinski gasket SGd^bin-) at stage n with dimension d equal to 
two, three and four for 6 = 2, and layer b equal to three for d = 2. The upper and lower bounds for 
the asymptotic growth constant, defined as zsCdb ~ li™i'->oo lii?TT-^^b(n)/f where v is the number 
of vertices, on these Sierpinski gaskets are derived in terms of the results at a certain stage. The 
numerical values of these zsg^ t evaluated with more than a hundred significant figures accurate. 
We also conjecture the upper and lower bounds for the asymptotic growth constant zsg^ 2 with 
general d. 
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I. INTRODUCTION 



The lattice gas with repulsive pair interaction is an important model in statistical me- 
chanics [l-4|. For the special case with hard-core nearest-neighbor exclusion such that each 
site can be occupied by at most one particle and no pair of adjacent sites can be simulta- 
neously occupied, the partition function of the lattice gas coincides with the independence 
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ynomial in combinatorics 



This model is a problem of interest in mathematics 



10|. While an activity (or fugacity) A can be associated to each occupied site, the spe- 
case with A = 1 counts the number of independent (vertex) sets Njs{G) on a graph 
Kaplansky considered the number of fc-element independent sets on the path and 



circuit graphs almost 70 years ago 



121 ] ■ For a graph G with v{G) vertices, the number of 



independent sets grows exponentially when v{G) is large. For the m xn grid graph, i.e. the 
square lattice (sq), it was shown that the limit \imrn,n->-oo Nis{sqy^™-^ exists and its upper 



and lower bounds were estimated 
been considered on various graphs 



13 [. The number of independent sets and its bounds had 
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16|. 



It is of interest to consider independent sets on self-similar fractal lattices which have 
scaling invariance rather than translational invariance 17|]. Fractals are geometric struc- 



tures of (generally noninteger) Hausdorff dimension realized by repeated construction of an 
elementary shape on progressively smaller length scales |l8, |l^. A well-known example 
of a fractal is the Sierpinski gasket which has been extensively studied in several contexts 
2o|-|36|. 

We shall derive the recursion relations for the numbers of independent sets on the Sier- 
pinski gasket with dimension equal to two, three and four, and determine the asymptotic 
growth constants. We shall also consider the number of independent sets on the generalized 
two-dimensional Sierpinski gasket with layer equal to three. 



II. PRELIMINARIES 



We first recall some relevant definitions for graphs and the Sierpinski gasket in this section. 
A connected graph (without loops) G = (V, E) is defined by its vertex (site) and edge (bond) 



sets V and E 



381]. Let v{G) = \ V\ be the number of vertices and e{G) = \E\ the number 



of edges in G. The degree or coordination number ki of a vertex Vi E V is the number of 
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edges attached to it. A /c- regular graph is a graph with the property that each of its vertices 
has the same degree k. An independent set is a subset of the vertices such that any two of 
them are not adjacent. 

When the number of independent sets Njs{G) grows exponentially with v{G) as v{G) — >■ 
oo, let us define a constant zq describing this exponential growth: 

Zg = lim 7777 — , (2.1) 

where G, when used as a subscript in this manner, implicitly refers to the thermodynamic 
limit. We will see that the limit in Eq. (12. ip exists for the Sierpinski gasket considered in 
this paper. 

The construction of the two-dimensional Sierpinski gasket SG2{n) at stage n is shown 
in Fig. [U At stage n = 0, it is an equilateral triangle; while stage (n + 1) is obtained by 
the juxtaposition of three ra-stage structures. In general, the Sierpinski gaskets SGd can be 



built in any Euclidean dimension d with fractal dimension D = ln{d + l)/ln2 2l|. For the 
Sierpinski gasket SGdin), the numbers of edges and vertices are given by 

v{SGd{n)) = ^[{d + ir + l] . (2.3) 

Except the {d + 1) outmost vertices which have degree d, all other vertices of SGd{n) have 
degree 2d. In the large n limit, SGd is 2(i-regular. 
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5G2(0) SG2{1) SG2{2) SG2{3) 
FIG. 1: The first four stages 71 = 0, 1, 2, 3 of the two-dimensional Sierpinski gasket SG2{n). 

The Sierpinski gasket can be generalized, denoted as SGd,b{n), by introducing the side 
length b which is an integer larger or equal to two [39]. The generalized Sierpinski gasket at 
stage n + 1 is constructed from b layers of stage n hypertetrahedrons. The two-dimensional 



SG2,b{n) with 6 = 3 at stage n = 1,2 and 6 = 4 at stage n = 1 are illustrated in Fig. |2j 
The ordinary Sierpinski gasket SGd{n) corresponds to the 6 = 2 case, where the index b is 
neglected for simplicity. The Hausdorff dimension for SG^.b is given by D = In / ln6 



39| . Notice that SGd,b is not A;-regular even in the thermodynamic limit. 




SG2,3{1) 



5^2,4(1) 



FIG. 2: The generalized two-dimensional Sierpinski gasket SG2fi{n) with 6 = 3 at stage n = 1, 2 and b = A 
at stage n = 1. 



III. THE NUMBER OF INDEPENDENT SETS ON 5^2 (n) 

In this section we derive the asymptotic growth constant for the number of indepen- 
dent sets on the two-dimensional Sierpinski gasket SG2{n) in detail. Let us start with the 
definitions of the quantities to be used. 

Definition III.l Consider the generalized two-dimensional Sierpinski gasket SG2,b{n) at 
stage n. (i) Define m2,b{n) = Njs{SG2,b{n)) as the number of independent sets, (ii) Define 
f2,b{n) as the number of independent sets such that all three outmost vertices are not in the 
vertex subset, (ii) Define g2,b{n) as the number of independent sets such that one certain 
outmost vertex, say the topmost vertex as illustrated in Fig. 0/or ordinary Sierpinski gasket, 
is in the vertex subset. (Hi) Define h2fi{n) as the number of independent sets such that two 
certain outmost vertices, say the leftmost and rightmost vertices as illustrated in Fig. [3 are 
in the vertex subset, (iv) Define P2,b{n) as the number of independent sets such that all three 
outmost vertices are in the vertex subset. 

Since we only consider the ordinary Sierpinski gasket in this section, we use the notations 
m2{n), f2{n), g2{n), h2{n) andp2(^) for simplicity. They are illustrated in Fig. |3l where only 
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the outmost vertices are shown. Because of rotational symmetry, there are three possible 
g2{n) and three possible h2{n) such that 

m2{Ti) = f2{n) + 35(2(72) + 3/i2(n) + p2{n) ■ (3.1) 

The initial values at stage zero are /2(0) = 1, 52(0) = 1, /i2(0) = 0, ^2(0) = and 777,2(0) = 4. 
The purpose of this section is to obtain the asymptotic behavior of m2{n) as follows. The 
four quantities f2{n), g2{n), h2{n) and P2{n) satisfy recursion relations. 

/!^\\ //^\k 

/2(n) 52 (n) /i2(n) P2{n) 

FIG. 3: Illustration for the configurations f2{n), g2in), h2{n), and P2{n). Only the three outmost vertices 
are shown explicitly, where a solid circle is in the vertex subset and an open circle is not. 

Lemma III.l For any non-negative integer n, 

f2{n + 1) = /|(n) + 3f2{n)gl{n) + 3gl{n)h2{n) + hl{n) , (3.2) 

g2{n + l) 

= flin)g2{n) + 2f2{n)g2{n)h2{n) + gl{n) + 2g2{n)hl{n) + 52(^)^2 (^) + hl{n)p2in\^ -3) 
h2{n + l) 

= f2{n)gl{n) + f2{n)hl{n) + 2gl{n)h2{n) + hl{n) + 2g2{n)h2{n)p2{n) + /72(n)p2(n)(5.4) 

P2{n + 1) = glin) + ?>g2{n)hl{n) + ?>hl{n)p2{n) + plin) . (3.5) 

Proof The Sierpinski gasket SG2{n + 1) is composed of three SG2{n) with three pairs 
of vertices identified. The number f2{n + 1) consists of (i) one configuration where all three 
SG2{n) belong to the class that is enumerated by f2{n); (ii) three configurations where one 
of the SG2{n) belongs to the class enumerated by f2{n) and the other two belong to the class 
enumerated by g2{n); (iii) three configurations where two of the SG2{n) belong to the class 
enumerated by g2{n) and the other one belongs to the class enumerated by h2{n)] (iv) one 
configuration where all three SG2{n) belongs to the class enumerated by h2{n) as illustrated 
in Fig. m Eq. (13.21) is verified by adding these configurations. 
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FIG. 4: Illustration for the expression of /2(n + 1). The multiplication of three on the right-hand-side 
corresponds to the three possible orientations of SG2{n + 1). 

Similarly, g2{n+l), h2{n+l) a.nd p2{n+l) for SG2{n+l) can be obtained with appropriate 
configurations of its three constituting SG2{n) as illustrated in Figs. [5], [6] and [7] to verify Eqs. 
(13. 3p . (13. 4p and (13. 5p . respectively. There are always 8 = 2'^ terms (counting multiplicity) in 
Eqs. (13. 2p - (13. 5p because for each of the three pairs of identified vertices it can be either in 
the vertex subset or not. □ 

^A^ A^ A'^A^A 

FIG. 5: Illustration for the expression of .92 (?i + !)■ The multiplication of two on the right-hand-side 
corresponds to the reflection symmetry with respect to the central vertical axis. 

^ A' A" A"" A' A"' A 

FIG. 6: Illustration for the expression of h2{n + 1). The multiplication of two on the right-hand-side 
corresponds to the reflection symmetry with respect to the central vertical axis. 

Alternatively, it is known that the number of dimer-monomers on a graph G is the same 
as the number of independent sets on the associated line graph L{G) [40]. Consider the 
sequence of graph H{n) shown in Fig. [H] that is obtained by adding an extra edge to each 
of the outmost vertices of the Hanoi graph. As H{n) has SG2{n) as its line graph, the 
enumeration of the number of independent sets on SG2{n) is equivalent to the enumeration 
of the number of dimer-monomers on these H{n). One can define corresponding quantities 
of f2{n), g2{n), h2{n), p2{n) on H{n) that satisfy the same recursion relations as in Lemma 
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FIG. 7: Illustration for the expression of P2{n + !)■ The multiplication of three on the right-hand-side 
corresponds to the three possible orientations of SG2{n + 1). 




HiO) H{1) 



FIG. 8: The first four stages n = 0, 1, 2, 3 of the graph H{n). 

The values of f2{n), (72(^)5 ^2(^)5 P2{n) for small n can be evaluated recursively by Eqs. 
(13. 2p - (13. 5p as listed in Table [H These numbers grow exponentially, and do not have simple 
integer factorizations. To estimate the value of the asymptotic growth constant defined in 
Eq. (12. ip . we need the following lemmas. 

For the generalized two-dimensional Sierpinski gasket SG2.b{n), define the ratios 



J2,b{'^) 92,b[n) h2,b{n) 

where n is a positive integer. For the ordinary Sierpinski gasket in this section, they are 
simplified to be a2{n), /32{n), 72 ('^) and their values for small n are listed in Table HTl From 
the initial values of f2{n), g2{n), h2{n), P2{n), it is easy to see that f2{n) > g2{n) > h2{n) > 
P2{n) for all non-negative n by induction. It follows that a2{n) , f32{n) , 'y2{n) G (0, 1]. 



Lemma III. 2 For any positive integer n, the ratios satisfy 

Oi2{n) < I32{n) < 72 (n) . 



(3.7) 



When n increases, the ratio a2{n) increases monotonically while 72 (?t-) decreases monotoni- 
cally. The three ratios in the large n limit are equal to each other 



lim a2{n) = lim /32(^) = hni 72 (^) • 



(3.8) 
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TABLE I: The first few values of f2in), g2{n), h2{n), P2{n)^ m2{n). 



n 





1 


2 


3 


4 


h{n) 


1 


4 


125 


4,007,754 


132,460,031,222,098,852,477 


92{n) 


1 


2 


65 


2,089,888 


69,073,020,285,472,159,669 


/i2(n) 





1 


34 


1,089,805 


36,019,032,212,213,865,476 


P2{n) 





1 


18 


568,301 


18,782,596,680,434,060,148 


m2{n) 


4 


14 


440 


14,115,134 


466,518,785,395,590,988,060 



TABLE II: The first few values of 02 (n), /32(n), 72(?^)• The last digits given are rounded off. 



n 


1 


2 


3 


4 


02 (n) 


0.5 


0.52 


0.521461147565444 


0.521463113425180 


Mn) 


0.5 


0.523076923076923 


0.521465743618797 


0.521463113431998 


l2{n) 


1 


0.529411764705882 


0.521470354788242 


0.521463113438816 



Proof It is clear that Eq. (13. 7p is valid for small values of n given in Table [TTl In order 
to save space, we will use «„, /3„, 7^ to denote a2{n), /32{n), 72(n) for the lengthy equations 
in this Lemma. By definition, we have 

Bn „ Cn. D' 



a„+i = , = a„— , 7„+i = a„— (3.9) 

■^n 



for a positive n, where 



An = 1 + 3a^ + 3a^/3„ + , 

Br, = 1 + 2anf3n + + l^lfil + O^Jnln + a^/3^7n , 

Dn = 1 + 3/31 + 3f3l^n + Phi , (3.10) 



such that 



ttn+i - ttn = -^|2a^(l + a„/3„)(/3„ - a„) + a^/3„(l + /3^)(7n - an)| , (3.11) 
/3n+i-an = -^|a„(/3n + an + an/3^ + an/5n7n)(/5n-an) + a^/3^(l + /3n7n)(7n-an)} • (3.12) 



It follows that 



/3n+l — Oin+1 



H ^ a 

In - 



2 



(3.13) 



where 



AA = 1 + + 2q;„^„ + 3a^ + 9q;^^„ + 2q;^/3^ + a^Mn + 3q;^^„ + Ua^^i + q;^/3^ 

+3a:t/3„7n + alPhn + 7a^^^ + 2^;^^^ + 3a^^^7n + ^o^^lln + 2a^^^ + 4a^^^7„ 

+«n/9n7n • (3-14) 

Using the fact that q;„, 7„ e (0, 1] and the inequality ^„ < 7^ to be shown below, q;„ < ^„ 
is proved by induction. Define = 7^ — «„, which is larger than 7^ — ^„ and ^„ — as we 
shall prove /3„ < 7n; then 

^2 . . 

^ r 1 



^2 



< ]4^{l + 3«n + 2an/9n + a'jl + + < 4 • (3-15) 
Similarly, we have 

7n+l - an = ^|2q;„/3„(1 + Pn'yn){Pn - Qin) + Qin/3^(1 + 7n) (7n - Qin) } , (3.16) 

and 

%--^|2q;„^„(1 + l3n1n)Bn - Otnil^n + + Oinl^l + Q;„/5„7„)C„| 

+ ^|^{«n/3^(l + 7^)i^n - a^/3^(l + Mn)Cn} 

- Oir. 



7n+l — Pn+1 



B„C, 



" " an/3nS (1 + OinPn)iPn " ^n) + ^n(7n " ^n/^n + 2Q;n/3n7n) (7n " ^n) 
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Wjnlnil + /3n7n)(7n " /3n)} , (3.17) 



where 



Using the fact that a„,/3„,7„ G (0, 1] and the inequahty a„ < /3„ given above, /3„ < 7„ is 
proved by induction. We also have 

7„+i - < — a„(l + Pnln + Pl + anPl + alPt) 

+a„/3^(l + an(3n + Pnln + 2a„/3.^7„ + a^/3„7„ + al(3^'jl)^ 
= «n + 2«„/3^ + «n/3„7n + 2a^/3^ + al(3^ + «„/3^7„ + «^/3^7„ 

■|l + 3a„/3„ + 3a„/3^ + 2alPl + 2a^/337„ + a3;34^^| < ^2 ^g^^g^ 



< 



BnCn 

From Eqs. fl3.15p and fl3.19p . we obtain e„+i < 2e^ for all positive n by induction. It follows 
that for any positive integer m < n, 

en < 2eti < 2[2el_,f <■■■< ^[2e^]'""". (3.20) 

Taking m as an integer larger than one, then the values of (3n, In are close to each other 
when n becomes large. 

Finally, it is clear that a2{n) increases monotonically as n increases by Eq. ( 13. lip . As 

7n-7n+i = ^{(l + /3n)(7n- an) + 2an/3„(l + /3„7„)(7„ -/?„)} , (3.21) 

we know 72 (n) decreases monotonically as n increases, and the proof is completed. □ 

Numerically, we find 

lim a2{n) = lim f32{n) = lim 72(n) = 0.521463113428094965776... (3.22) 

n— >oo n— >cxD n— >oo 

From the above lemma, we have the following bounds for the asymptotic growth constant. 
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Lemma III. 3 The asymptotic growth constant for the number of independent sets on 
SG2{n) is hounded: 

^ ln[/2(m)] + ^ ln[l + al{m)] < zsg, < ^ Hh{m)] + ^ ln[l + ^l{m)] , (3.23) 
where m is a positive integer. 

Proof From Eq. fl3.2p and Lemma [III. 2[ we have the upper bound for f2{n), 

f^(n) = fl{n - 1)[1 + 3a^(n - 1) + 3a^(n - l)^2(n - 1) + al{n - l)(3l{n - 1)] 

< fl{n - + ^l{n - l)f 

< {f^{n - 2)[1 + ^lin - 2)]^}'[1 + ^i{n - 1)]' < ■ ■ • 

< [f2im)r-[l + jlim)f'—''K (3.24) 
From Eq. (13. ip . the number of independent sets has the upper bound 



'm2{n) = f2{n)[l + 3a2{n) + 3a2{n)/32{n) + a2{n) I32{n)^2{n)] 

< [f2{m)r-[l + ^l{m)]l^'---%+^2{n)f . (3.25) 

As the number of vertices of SG2{n) is 3(3" + l)/2 by Eq. (12.31) . the upper bound for 2:5^2 
defined in Eq. (12. ip follows. The lower bound for zsg2 can be derived similarly. □ 

As m increases, the difference between the upper and lower bounds in Eq. (I3.23P becomes 
small and the convergence is rapid. The numerical value of zsg2 can be obtained with more 
than a hundred significant figures accurate when m is equal to eight. 

Proposition III.l The asymptotic growth constant for the number of independent sets 
on the two-dimensional Sierpinski gasket SG2{n) in the large n limit is zsg2 = 
0.38430953443368558352.... 

As mentioned previously, the number of dimer-monomers on the graph H{n) illustrated 
in Fig. [8] is the same as the number of independent sets on the two-dimensional Sierpinski 
gasket SG2{n). Similar to Eq. (12. ip . one can define a constant for the exponential growth 
of the number of dimer-monomers: 

l-'^^^f ) . (3.26) 

where N£,m{G) is the number of dimer-monomers on a graph G. As the number of vertices 
of H{n) is 3" + 3, we have the following corollary. 
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Corollary III.l The asymptotic growth constant for the number of dimer-monomers on 
the graph H{n) in the large n limit is z'^ = 0.57646430165052837528.... 

To the best of our knowledge, this result has not been obtained elsewhere. 

IV. THE NUMBER OF INDEPENDENT SETS ON 6^^2,3 (n) 

The method given in the previous section can be applied to the number of independent 
sets on SGd,b{n) with larger values of d and b. The number of configurations to be considered 
increases as d and b increase, and the recursion relations must be derived individually for 
each d and b. In this section, we consider the generalized two-dimensional Sierpinski gasket 
SG2,b{n) with the number of layers b equal to three. For SG2,zij'), the numbers of edges 
and vertices are given by 

e{SG2M)) = 3 X 6^ (4.1) 
7 X 6" -I- 8 

v{SG2M) = -—- , (4.2) 

5 

where the three outmost vertices have degree two. There are (6" — l)/5 vertices of >S'G'2,3('^) 
with degree six and 6(6" — l)/5 vertices with degree four. The initial values for the number 
of independent sets with various conditions are the same as those for SG2'- f 2,3(0) = 1, 
^72,3(0) = 1, /i2,3(0) = and ^2,3(0) = 0. The recursion relations for SG2,3{n) are lengthy and 
given in the appendix. Some values of f2,3{n), g2,3{n), h2,^{n), p2,3{n), m2,z{n) are hsted in 
Table IIIII These numbers grow exponentially, and do not have simple integer factorizations. 



TABLE III: The first few values of /2,3(n), g2,3{'n), h2,3{n), P2,3in), m2,3{n). 



n 





1 


2 


3 


f2,3{n) 


1 


19 


172,371,175 


93,818,345,014,803,648,739,612,995,034,820,933,103,277,876,214,071 


52,3 (n) 


1 


9 


80,291,169 


43,700,938,182,461,202,772,695,141,988,444,331,720,442,482,282,619 


h2,3{n) 





4 


37,399,906 


20,356,061,468,851,869,739,344,457,713,631,919,274,541,443,648,604 


P2,3(ri) 





2 


17,420,990 


9,481,930,039,890,479,716,613,035,420,873,292,623,048,215,623,126 


m2,3(ra) 


4 


60 


542,865,390 


295,471,274,008,633,345,992,344,829,561,922,978,711,277,869,630,866 
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The values of the ratios 02,3 ('^), f32,3{n), 72,3 (^) defined in Eq. (13. 6p for small n are listed 
in Table ITVl The sequence of 0:2,3 (^) decreases monotonically as n increases, while P2,3{n) 
increases monotonically. Except the first term 72,3(1), 72,3('^) also increases monotonically 
for n>2. We again have 02,3 1?^), /32,3('^), 72,3 (^) ^ (0, 1] but 72,3(^) < P2,3{n) < a2,3{n) for 
n > 2, in contrast to Lemma [III. 2[ 

TABLE IV: The first few values of 0:2,3 (?^), I32,3in'), 72,3 (?i-)- The last digits given are rounded off. 



n 


1 


2 


3 


a2,3(«) 


0.47368421052631578947 


0.46580391994195085112 


0.46580376338514186621 


^2,3 (n) 




0.46580348082863259844 


0.46580376338514186620 


72,3 (?^) 


0.5 


0.46580304239267339335 


0.46580376338514186619 



By the same argument given in Lemma [111.31 we have the upper and lower bounds of the 
asymptotic growth constant for the number of independent sets on SG2,3{n): 

,^^|51n/2,3(m) + In [1 + jl,{m)] + 61n [1 + 72,3('^)]} < zsg,,s 

- 7^{51n/2,3M + In [1 + al,im)] + 61n [1 + al,{m)]] , (4.3) 

with m a positive integer. The convergence of the upper and lower bounds remains rapid. 
More than a hundred significant figures for zsg2 3 can be obtained when m is equal to five. 
We have the following proposition. 

Proposition IV. 1 The asymptotic growth constant for the number of independent sets 
on the two-dimensional Sierpinski gasket SG2,3{n) in the large n limit is zsg2,3 — 
0.38135033366164857274.... 

V. THE NUMBER OF INDEPENDENT SETS ON SG3{n) 

In this section, we derive the asymptotic growth constant of independent sets on the 
three-dimensional Sierpinski gasket SG3{n). We use the following definitions. 
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Definition V.l Consider the three-dimensional Sierpinski gasket SG^ln) at stage n. (i) 
Define m^{n) = Nis{SG^{n)) as the number of independent sets, (ii) Define fsin) as the 
number of independent sets such that all four outmost vertices are not in the vertex subset. 
(Hi) Define g^{n) as the number of independent sets such that one certain outmost vertex is in 
the vertex subset, (iv) Define h^{n) as the number of independent sets such that two certain 
outmost vertices are in the vertex subset, (v) Define p^{n) as the number of independent 
sets such that three certain outmost vertices are in the vertex subset, (vi) Define qs{n) as 
the number of independent sets such that all four outmost vertices are in the vertex subset. 



The quantities fsi^n), g^iji), h^{n), p^in) and q^in) are illustrated in Fig. [9l where only 
the outmost vertices are shown. There are (^) = 4 equivalent gsin), (2) = 6 equivalent 
h3,{n), and (^) = 4 equivalent Pz{n). By definition, 

m^{n) = fsin) + igs^n) + Gh^i^n) + 4p3(n) + q^^n) . (5.1) 

The initial values at stage zero are /3(0) = 1, ^'3(0) = 1, /i3(0) = 0, ^3(0) = 0, ^3(0) = and 
^3(0) = 5. 








fsin) gsin) hz{n) pz{n) qz{n) 

FIG. 9: Illustration for the configurations fs^n), g^in), h^^n), p^in) and q^in). Only the four outmost 
vertices are shown explicitly, where a solid circle is in the vertex subset and an open circle is not. 

The recursion relations are lengthy and given in the appendix. Some values of fsin), 
g^in)^ h^^n), Pziji), qs^n), m^{n) are listed in Table El These numbers grow exponentially, 
and do not have simple integer factorizations. 

Define ratios 

( \ 93{n) . . . hs{n) p^{n) q^{n) 

/3(^) 93{n) h^[n) p^[n) 

for a positive integer n as in Eq. f l3.6p . As n increases, we find a^{n) increases monotonically 

while Pziji), 73 (^), '^3('^) decreases monotonically with the relation 03 (n) < Pziji) < 73 (^) < 

5z{n). The values of these ratios for small n are listed in Table IVTl Numerically, we find 

lim a3(n) = lim (3^{n) = lim 73(n) = lim S^in) = 0.442256573677178603386... (5.3) 
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TABLE V: The first few values of fz{n), gz{n), h^{n), pz{n), qz{n), msin). 



n 





1 


2 


3 


h{n) 


1 


10 


25,817 


1,292,964,293,737,151,090 




1 


4 


11,387 


571,820,791,550,665,532 


hsin) 





2 


5,050 


252,892,039,471,313,074 







1 


2,252 


111,843,868,747,687,217 







1 


1,010 


49,464,202,269,253,193 


m3(n) 


5 


43 


111,683 


5,594,439,374,027,693,723 



TABLE VI: The first few values of 03(71-), /S^^n), 73(71-), 5s{n). The last digits given are rounded 
off. 



n 


1 


2 


3 


4 


03(71) 


0.4 


0.441065964287098 


0.442255671189410 


0.442256573676665 




0.5 


0.443488188284886 


0.442257510059261 


0.442256573677711 


l3{n) 


0.5 


0.445940594059406 


0.442259349015113 


0.442256573678758 




1 


0.448490230905861 


0.442261188057088 


0.442256573679804 



By a similar argument as Lemma lin.31 the asymptotic growth constant for the number 
of independent sets on SGsi^n) is bounded: 

^ HMm)] + ^ ln[l + alim)] < zsg. < ^ ln[/3(m)] + 1, ln[l + 6l{m)] , (5.4) 

where m is a positive integer. More than a hundred significant figures for zsgs can be 
obtained when ttt, is equal to seven. We have the following proposition. 

Proposition V.l The asymptotic growth constant for the number of independent sets 
on the three-dimensional Sierpinski gasket SG^^n) in the large n limit is zsgs = 
0.32859960572147955761.... 
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VI. THE NUMBER OF INDEPENDENT SETS ON SGi{n) 



In this section, we derive the asymptotic growth constant of independent sets on the 
four-dimensional Sierpinski gasket SGi^{n). We use the following definitions. 

Definition VI. 1 Consider the four-dimensional Sierpinski gasket SGi{n) at stage n. (i) 
Define m4^{n) = Njs{SG4^{n)) as the number of independent sets, (ii) Define fi{n) as the 
number of independent sets such that all five outmost vertices are not in the vertex subset. 
(Hi) Define gi{n) as the number of independent sets such that one certain outmost vertex is in 
the vertex subset, (iv) Define h^ln) as the number of independent sets such that two certain 
outmost vertices are in the vertex subset, (v) Define p4{n) as the number of independent sets 
such that three certain outmost vertices are in the vertex subset, (vi) Define qi{n) as the 
number of independent sets such that four certain outmost vertices are in the vertex subset, 
(vi) Define ri{n) as the number of independent sets such that all five outmost vertices are 
in the vertex subset. 

The quantities fiiji), g4{n), h^^iji), p^i^n), qiin) and r4(n) are illustrated in Fig. [101 where 
only the outmost vertices are shown. There are (^) = 5 equivalent (n), (2) = 10 equivalent 
/i4(n), (3) = 10 equivalent P3(n), and (^) = 5 equivalent q^in). By definition, 

m^^n) = fiin) + ^g^in) + 10/14(72) + IQp^in) + 5g4(n) + r^s^in) . (6.1) 



The initial values at stage zero are /4(0) = 1, g^ifS) = 1, /i4(0) = 0, ^4(0) = 0, ^4(0) = 0, 
r4(0) = and 7714(0) = 6. 




fiin) gi{n) h/^{n) pi{n) 94(71) r4(7i) 

FIG. 10: Illustration for the configurations f^in), gi{n), hi^n), P4(n), (74(11) and r4(ii). Only the five 
outmost vertices are shown explicitly, where a solid circle is in the vertex subset and an open circle is not. 



The recursion relations are lengthy and given in the appendix. Some values of fiin), gi{n), 
hiin), P4{n), q4{n), r^i^n), m^^n) are listed in Table lylTI These numbers grow exponentially, 
and do not have simple integer factorizations. 
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TABLE VII: The first few values of fi{n), g4{n), hi{n)^ Piin), 54 (n), r4(n), 7714(71). 



77 





1 


2 


3 


Uin) 


1 


26 


48,645,865 


1,209,689,823,065,753,613,801,849,265,389,348,210,254 


54 N 


1 


10 


19,499,025 


485,275,031,765,121,996,003,377,748,244,728,141,942 


h/i{n) 





4 


7,827,058 


194,671,321,306,020,419,533,199,834,929,606,628,798 


Pi{n) 





2 


3,146,558 


78,093,721,039,746,646,163,976,217,053,630,607,240 


q4{n) 





1 


1,266,948 


31,327,833,873,772,900,771,790,623,812,192,536,505 







1 


510,980 


12,567,379,442,065,248,794,102,222,711,306,394,841 


7774(77) 


6 


142 


262,722,870 


6,532,921,954,159,964,003,443,553,868,217,630,357,710 



Define ratios 

g^in) h^in) p^in) q^in) n{n) 

(6.2) 

for a positive integer n as in Eq. ( 13.60 . As 77, > 2 increases, we find 04 (n) increases 
monotonically wliile ^^{n), 74 (^), S^ln), r]4{n) decreases monotonically with the relation 
o:4^{n) < (34^{n) < 74 (n) < ^4(77) < rji{n). The values of these ratios for small n are listed in 
Table IVIIIi Numerically, we find 

lim aA{n) = hm /S^in) = lim 74(77) = lim 64(71) = lim r]4in) = 0.401156636030339443965. 

(6.3) 

By a similar argument as Lemma IIII.31 the asymptotic growth constant for the number 
of independent sets on SGii^n) is bounded: 

^ ln[/4(m)] + ^ ln[l + al{m)] < zsg, < ^ HU^)] + ^ + vlim)] , (6.4) 

where m is a positive integer. More than a hundred significant figures for zsd can be 
obtained when 777 is equal to seven. We have the following proposition. 



Proposition VI. 1 The asymptotic growth constant for the number of independent sets 
on the four- dimensional Sierpinski gasket SG^ln) in the large n limit is = 
0.28916553234872775551.... 
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TABLE VIII: The first few values of 04 (n), /34(n), 74 (n), 5i{n), ?/4(n). The last digits given are 
rounded off. 



n 


1 
i 



L 


Q 


4 


04(71) 


0.384615384615385 


0.400836227292906 


0.401156579572832 


0.401156636030338 


/34(n) 


0.4 


0.401407660126596 


0.401156681393497 


0.401156636030341 


74 (n) 


0.5 


0.402010308343186 


0.401156783217105 


0.401156636030344 




0.5 


0.402645684586141 


0.401156885043655 


0.401156636030347 


rj4{n) 


1 


0.403315684621626 


0.401156986873147 


0.401156636030350 



VII. BOUNDS OF THE ASYMPTOTIC GROWTH CONSTANTS 

For the rf-dimensional Sierpinski gasket SGd{n), we conjecture that similar upper and 
lower bounds for the asymptotic growth constant as in Lemma IIII.3I hold, 

2 1 

Hfd{m)] + jj—^\n[l + «^(m)] < zsg, 

< (^TT)^ HhM] + ^ ln[l + dm)] (7.1) 

with a positive integer m, where the ratios are defined as 

/ N 9d{n) td{n) 
ad{,n) = —— , Cd{n) = —— , 7.2 
fd{n) Sd{n) 

for a positive integer n. fd{n) again is the number of independent sets such that all d + 1 
outmost vertices are not in the vertex subset, Qdin) is the number of independent sets such 
that one certain outmost vertex is in the vertex subset, Sd{n) is the number of independent 
sets such that all but one certain outmost vertex are in the vertex subset, and t^ln) is the 
number of independent sets such that all c? + 1 outmost vertices are in the vertex subset. 

Although the quantities in Eq. (17. ip for general m is difficult to obtain, one can consider 
the simplest case m = 1. Denote the upper and lowers bounds at m = 1 as zsg^ ^SGdi 
respectively. Because 5^(1) = td(l) = 1 and Qdi)-) = fd~i{^), we have 
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TABLE IX: Numerical values oi z_gQ^, ^SGa^ some ratios of them to zsca- The last digits given 
are rounded off. 



d 




l5Gd 




ZSGi 




^ScJzsGi 


2 


4 


0.3824465974 


0.5391144738 


0.3843095344 


0.9951525088 


1.402813164 


3 


10 


0.3249281379 


0.4611099318 


0.3285996057 


0.9888269257 


1.403257715 


4 


26 


0.2882396119 


0.3992771592 


0.2891655323 


0.9967979570 


1.380790981 


5 


76 


0.2590427565 


0.3561208268 








6 


232 


0.2368781125 


0.3213368369 








7 


764 


0.2184809121 


0.2940986410 








8 


2620 


0.2034116955 


0.2713602941 








9 


9496 


0.1905090814 


0.2524872368 








10 


35696 


0.1794854089 


0.2362827010 









and the task reduces to the determination of fd{^)- It is easy to see that /i(l) = 2 and 
assume the formal quantity /o(l) = 1, then /d(l) satisfies the recursion relation 

/,(!) = /,_i(l) + rf/,_2(l) (7.4) 

for d > 2. This relation can be understood as follows. The d- dimensional Sierpinski gasket 
SGd{l) at stage one is the juxtaposition of d+1 complete graphs -ft'rf+i. For the enumeration 
of /d(l), consider one of the complete graphs. In the case that all d interior vertices of the 
complete graph are not in the vertex subset, the number is the same as Qdi^) = /d-i(l), 
which is given as the first term on the right-hand-side of Eq. (17.40 . In the case that one 
of the d interior vertices of the complete graph is in the vertex subset, the number is given 
by /d-2(l), which gives the second term on the right-hand-side of Eq. (17.41) . It follows 
that /d(l) is equal to the number of permutation involutions on d+1 elements. The values 
of /d(l), ZgQ^, zsGd for small d are listed in Table [IXl We notice that z^q^ is closer to 
zsGd compared with zsg^ and serves as an approximation for zsg^- The fact that zsgj, 
decreases as d increases is easy to understand because the degrees of the vertices of SGdiji) 
also increase. 
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Appendix A: Recursion relations for SG2,3{n) 



We give the recursion relations for the generalized two-dimensional Sierpinski gasket 
SG2,3{n) here. Since the subscript is d = 2, 6 = 3 for all the quantities throughout this 
section, we will use the simplified notation to denote f2,3{n + 1) and similar notations 
for other quantities. For any non-negative integer n, we have 



,27,4 I 7,6 I t3„3 I cJ-2„3i, , at „3i,2 , o-f2„ 7,3 , o -f „4„ , o„3i,3 



+^fn9lhlpn + QgnKPn + Qfn9nhl + Sfnh^Pn + 9c/^/i^Pn + ^9nhnPl + ^9nhlpl 

+hlpl , (Al) 

gn+i = fn9n + 2/,t^„/i„ + Af^g^ + ifngl + ^ flglK + 2/^^„/i^ + flglpn + AglK 

+10fnglhl + 2flgnhl + 2/„^>„ + 2/^^^/i„p„ + Iglhl + 3/„^^/i^p„ + Sg^KPn 
+2UgnK + '^9nhl + Aglhlp^ + + flgt, + 4/„^>„ + 2flglhl + 3^^/i^ 

+8fn9lhl + 2fr,glKPn + fhuh^Pn + ^^n + 4y>^ + Afn9nhlPn + 8^^/l^„ 

+2fnglKpl + 2fnhl + Ignhipn + 5^>^p^ + 2/„/i^^ + glpl + 2/1^^ ^ 4^^/,2^3 

, (A2) 

= fn9l + 2/'(7' + AflglK + 7/„(7^/i„ + 5/^(7^/1^ + 2flglpn + /^^ ^ 4^4^2 
+8fnglhl + AfnglKpn + 2f^gnhlpn + 2glpn + 3^^/i^ + 2/„^„/i^p„ + Sglhlpn 

+7^^^^ ^ Sf^glhlp^ + 2g^hnPn + ^fn9nK + 45n/in + + 10^>^„ 

+2^n^nPn + '^fn9nhlpl + 9y„/i^^ + 3/i^„ + 2glKpl + /n/i^^ + 4/i^^ + 2g^Kpi 
+hnPl , (A3) 

= + gy2^3^^ ^ 9^^^3^2 ^ 3^2^^^3 + ^f^^^ + 253^3 + e/,^^^^^ ^ g^4^^^^ 
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+9^>^ + Qglhlp^ + 2hl + l2gXnPn + ^glhlpl + 9/i:^p2 ^ Qg^h^pl + 6/i^^ + . 

(A4) 

There are always 128 = 2^ terms (counting multiplicity) in these equations. 
Appendix B: Recursion relations for SG2,{n) 

We give the recursion relations for the three-dimensional Sierpinski gasket SG^in) here. 
Since the subscript is c? = 3 for all the quantities throughout this section, we will use the 
simplified notation f^+i to denote f^{n + 1) and similar notations for other quantities. For 
any non-negative integer n, we have 

= + Qflgl + UUglK + K + ^fnhl + I2glhl + 45^„ + ?,hi + 125„/iX 
+6/iX + Pn, (Bl) 

5'n+l = flgn + ^fngnhn + ^Ugl + Qfngnhl + Gg^hn + ^fnglPn + T gnh^ + Sfnhlpn 

+9glhnPn + glqn + 6/l^„ + QgnKpl + ^gnKqn + ^KpI + '^hlpuQu + P^Qn (B2) 

K+l = fWn + ^fnglK + /n/in + + Vnhl + 7^^/l^ + 2^5?„ + Af^gnKPn + 2/l^ 

+12fl'n^nPn + 2glpl + 2fnhnpl + 2glhnqn + ^KpI + 25(„p^ + 2hlqn + 4:gnhnPnqn 

+Pn + 4:hnPlqn + /^n?^ + pIqI > (B3) 

Pn+1 = fng'i + ^fngnhl + ^glK + Qgnhi + Sfnhlpn + QglKPn + 7/l^p„ + '^gnKpl + /nPn 
+^gnhlqn + QhnPl + Ghlp^qn + ^gnPlqn + ^P^qn + "^KPriql + Pn?^ , (B4) 

+6p^?^ + q'n ■ (B5) 
There are always 64 = 2^ terms (counting multiplicity) in these equations. 

Appendix C: Recursion relations for SG4{n) 

We give the recursion relations for the four-dimensional Sierpinski gasket SG4{ri) here. 
Since the subscript is d = 4 for all the quantities throughout this section, we will use the 
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simplified notation fn+i to denote /4(n + 1) and similar notations for other quantities. For 
any non-negative integer n, we have 



+10p2g3 + ^5^ (CI) 



^n+i = fn9n + 4/„3^„/i„ + 6flgl + 12/„2^„/i^ + 24/„^^/i„ + Qflg^pn + 3^^ + 28/„^„/i^ 

+6fX^Pn + SGfnglhnPn + ^fnQnQn + SGg^hl + 10(/>„ + 24fnhlpn + 24/„^„/l„p^ 

+12fngnhlqn + Slfif^/i^ + OOg^hlpn + l^glpl + IQglKqn + S-^Vn + ^^frXiPl 
+12fnhlpnqn + 36/?.f,p„ + 102gnhlpl + 125(^pf, + 'i^glKpuQu + 3Qgnhiqn + Qg^hlrn 
+^fnPl(ln + + 52hlpnqn + GOgrAiPlQu + 3/?.f,r„ + Ug^hlpnrn + l^gnPn 

+12gnhlql + 52hnplqn + ^.Ohlp^ql + 12/iX^n + + 4/i3g„r„ + 12^„p^g^ 

+4c?nPn'^n + IS/lnPn^n + 12hnplqnrn + 18p^?n + SPn'^n + 4p„?^ + Gp^^^r^ + q^Tn , 

(C2) 



^n+l = /nfn + ^IIqIK + /n^n + 3/n5n + 3/,^/?;^ + 21 fng^hl + 6/„^j^p„ + QflgnKPn 
+^gtiK + + Sdfngnhlpn + ^fnOlpl + '^flKvl + 31^^/i^ + QfuQlKQu 

+ZOglKPn + 2^^?n + 21/n/lX + + 6/„/l^?n + 12/n^nMn?n + 7/l^ 

+72gnhlpn + hlglKpl + 27 glhlqn + QglpuQu + 2g^Krn + Sfnpt, + 12fnhnplqn 

+^fnhlqn + 5Ahlpl + bAg^Kpl + 15/i^g„ + TSg^hlpnqn + QglKql + l^glplqn 
+6gnhlr,, + QglKPnTn + ^fnvWn + 27/;,„pf, + Slhlplqn + SOgrAiPnql + 13/i^g,^, 

+2QgnPUn + 12hlpnrn + 12gnKplr\, + GgnhlquTn + ^IhnPWn + 18Pn?n + 
+^9nPnqn + + IS/i^Png^r + GgnPlquTn + ^Kqn + ISp^?^ 

+ 12hnPnqlrn + ^KpItI + Up^q^Tn + + 6pn?n'^n + 3p^?nr^ (C3) 



Pn+l = /^^j^ + ^fndnhl + QfnglK + ^f, + 12/„^„/l^ + Sf^hlpn + 12 fnglKpn + 

+3fi'^Pn + Ufnh^Pn + ISfngnKpl + /nPn + '^fngnKqu + IS^-n/l^ + blg^hlpn 
+%5^n + ^glhnqn + 12fnhnPl + 12 f^hlp^qn + Qfn9nPlqn + 27/i^„ + Slg^hlpl 

+13^5^^ + SOglKpnqn + 26^„/i^gn + Sglhlvn + GfnPlqn + QfnKPnql + 54.hlpl 
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+27gnplql + ^dnKql + 125(„/i„p„g„r„ + GgnP^Vn + SOhnPuql + SGhnplquTn 

+31p^g^ + QhlqlVn + QPnVn + QgnPuqlVn + 2^„g,t + ^hlpnVl + ^KqlVn + Opn^n 

+21p^g^r„ + GhnPnqnrl + 3p^r^ + 3g^r„ + 6p„g^r^ + p^r^ + g^r^ , (C4) 

+4^fn9nPl + ^fnhlqn + 6/i^ + 52^„/i^p„ + 30^^/i„p^ + 12^>^g„ + Zfnpi 

+12fnhnPlqn + 54/?,^^^ + 52gnhnPl + ISh^qn + GOgnhlpnqn + I2glplqn + ^nhlr^ 

+^fnPlql + "i^KPn + 102/l^p2g^ _^ SQgnhnPnql + + ^QgnPlqn + I2hlpnrn 

+12g^hnPlrn + fnqn + QOhnPlql + 31p^g„ + Uhlqf^ + IQgnPnql + 2AKplrn 
+24:hlpnqnrn + I2gnplqnrn + lO/in^n + 36p^g^ + 36hnPnqlrn + Ghnplrl + 28plqnrn 
+4^n?n^n + 3g^ + 6/i„g^r^ + 24p„g^r„ + Uplq^rl + 6g^r^ + 4p„?„r^ + g^r^^ , (C5) 

Tn+i = 9n + ^Oglhl + I5gnh'^ + SOg^hlpn + iOh'^Pn + QOgnhlpl + lOfif^p^ + 20gnhlqn 
+70hlpl + GOhlpnqn + Q^gnKpiqu + 5/i^r„ + IhgnP^ + I20hnplqn + 60/i^p„g^ 
+30hlplrn + I2pl + SOgnPlql + QOKpuql + QOKplquTn + 70plql + ^^p^rn 
+5gnqt, + 2QKqlrn + 30p„?^ + 60p^g^r„ + lOp^r^ + 15g^r„ + 30p„g^r^ + lOg^r^ 
■ (C6) 

There are always 1024 = 2^° terms (counting multiplicity) in these equations. 
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